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Abstract. Mandelbrot introduced the concept of fractals to describe the non- 
Euclidean shape of many aspects of the natural world. In the time series context 
he proposed the use of fractional Brownian motion (fBm) to model non-negligible 
temporal persistence, the "Joseph Effect"; and Levy flights to quantify large dis- 
continuities, the "Noah Effect". In space physics, both effects are manifested in the 
intermittency and long-range correlation which are by now well-established features 
of geomagnetic indices and their solar wind drivers. In order to capture and quantify 
the Noah and Joseph effects in one compact model we propose the application of 
the "bridging" fractional Levy motion (fLm) to space physics. We perform an initial 
evaluation of some previous scaling results in this paradigm, and show how fLm can 
model the previously observed exponents. We suggest some new directions for the 
future. 

Keywords: 

1. Introduction 

Ever since it became clear that Earth's magnetosphere is influenced 
by the sun, significant effort has been devoted to establishing the re- 
lationship between fluctuations in the energy delivered by the solar 
wind to the magnetosphere and variations in the magnetospheric re- 
sponse. A particularly important diagnostic for the response has been 
the family of geomagnetic indices, especially the Auroral Electrojet 
index AE (Davis and Sugiura, 1966). A common proxy for the solar 
wind input is the e function (Perreault and Akasofu, 1978) which esti- 
mates the fraction of the solar wind Poynting flux through the dayside 
magnetosphere. 

One approach is to investigate causal relationships, and consider- 
able sophistication has now been developed in this (e.g. Ukhorskiy 
et al.,2004;March et al,2005 and references therein). However, even 
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without examining causality, significant information can be obtained 
by examining the scaling behaviour of fluctuations. A first analysis 
of this, in the Fourier domain, was done by Tsurutani et al.(1990) 
using the power spectrum. Subsequent analyses have introduced other 
methods for detecting scale invariance (e.g. Takalo et al.,1993;Freeman 
et al., 2000a, 2000b;Uritsky et al.,2001). Most recently Hnat et al., 
(2002a; 2002b; 2003a; 2003b; 2005) and Chapman et al.,(2005) have 
studied the scaling collapse of the increments of time series. 

A fundamental problem has been raised by the evidence for multi- 
fractality in some solar wind quantities (e. g. Hnat et al., (2002a) and 
references therein) and the AE index (Consolini et al., 1996). Multifrac- 
tality is physically well motivated-at least for solar wind quantities-in 
that it arises naturally from the intermittency of multiplicative turbu- 
lent cascade models (Frisch, 1995). Multifractality would imply that the 
Hurst's "roughness" exponent H is not constant but varies from scale 
to scale. This evidence for multifractality in the indices thus means that 
any comparison of pairs of scaling exponents derived from solar wind 
and geomagnetic indices may be problematic (Watkins, 2002; Chang 
and Consolini, 2001). Preliminary comparisons of solar wind and geo- 
magnetic field measurements made using multiscaling measures (Voros 
et al., 1998) showed similarity at low orders after low pass filtering 
of the magnetospheric quantities. However, Hnat et al. (2002-2005), in 
examining a range of solar wind quantities, have recently found some 
apparent simplifications. They see the intriguing result that although 
some quantities (notably v and B) do not show a simple scaling collapse, 
consistent with their well-known multifractality, others (such as B 2 ) 
do i.e. they are, in this sense, effectively monofractal. Recently Hnat 
et al. (2003b) have extended the 1 year AE/U/L dataset studied by 
Hnat et al.(2002b)to the 10 years used by Freeman et al. (2000a). They 
find that when such long auroral index datasets are examined, AE and 
e do indeed have discernably different PDFs. 

Such analyses are not easy to compare. Some used overlapping in- 
dex and solar wind time series (Uritsky et al., 2001), other did not 
(Freeman et al., 2000a). Techniques which impose finite limits on the 
integral used to evaluate structure functions have also been explored 
((Chapman et al., 2005) and references therein). The choice of solar 
wind measures and geomagnetic time series has also varied. It seems 
to us thus imperative to try to start to reconcile the various studies 
and understand why some show much greater similarity between the 
solar wind signal and indices than others. We also believe that the 
synthesis of observations will help towards a goal we have proposed else- 
where: The definition of models which are either I) simple, statistical, 
"strawman" models which may nonetheless capture some relevant fluc- 
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tuation phenomenology e.g. the fractional lognormal model sketched 
by Watkins(2002)) or II) more clearly statistical physics-based e.g. the 
generalised Fokker-Planck model of Hnat et al.(2005) and Chapman 
et al.(2005)). 

By analogy with mathematical economics we may think of the Type 
I models as modelling the "stylized facts" of the coupled solar wind 
magnetospheric system (Watkins, 2002). In this paper we shall in- 
troduce one such model: fractional Levy motion (Mandelbrot, 1995; 
Chechkin and Gonchar, 2000b), in order to see how well it can de- 
scribe the solar wind e function and the AE family of indices (AE 
itself, AU and AL). Preliminary comparison is made with some of the 
measurements listed above, and it is shown that the model provides 
a good quantitative explanation for the difference between two scaling 
exponents first noted in this context by Hnat et al. (2002a) as well as a 
possible qualitative explanation for the multifractal behaviour seen by 
Hnat et al. (2003b). Where relevant, the effect of the truncation (finite 
variance) implicit in both a natural data series and a computer model 
are noted. Future directions are then sketched. 



2. Datasets used 

The AE and e data are a 1 year subset of those studied by Hnat 
et al. (2002b; 2003a). They correspond to the years 1978 and 1995 re- 
spectively. As with (Hnat et al., 2002b) solar wind data is taken only for 
periods when WIND is definitely in the solar wind, see (Freeman et al., 
2000a) for details. We follow Hnat et al. (2002b) by firstly differencing 
the time series X(t) of the indices AE, AU, AL and e at intervals r 
of 1,2,3. . . times the fundamental sampling period (1 minute for the 
indices and and 46 seconds for e) to generate difference time series 
5X(t,r) = X(t + r) — X(t). For further details of the dataset and 
preprocessing techniques see (Hnat et al., 2002b) and references therein. 

3. Motivation for and testing of a fractional Levy motion 

model 

3.1. Fractional Levy motion as a bridge between Levy 
flights and fractional brownian motion 

As noted by Mandelbrot(1995): 

The "normal" model of natural fluctuations is the Wiener Brownian 
motion process (WBm). By this standard, however, many natural 
fluctuations exhibit clear-cut "anomalies" which may be due to 
large discontinuities ("Noah Effect") and/or non-negligible global 
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statistical dependence ("Joseph Effect"). [Mandelbrot's book "The 
Fractal Geometry of Nature"] ... shows that one can model various 
instances of the Noah effect by the classical process of [standard 
Levy motion] (sLm), and various instances of the Joseph effect by 
the process of [fractional Brownian motion] (fBm). 

Takalo et al.(1993) were the first to use fBm as a model of the auroral 
indices, but it subsequently could not describe the highly non-Gaussian 
leptokurtic distributions seen in differenced solar wind and geomagnetic 
index quantities. This can for example be seen in Fig. 7 of (Chapman 
et al., 2005) where the pdf of differences 5X of AE is contrasted with 
the Gaussian pdf of an fBm with equal Hurst exponent H. Similarly 
we are are aware of only a small number (Kabin and Papitashvili, 
1998; Consolini et al., 1997; Hnat et al., 2002a; Bruno et al., 2004) of 
discussions of the use of truncated sLm as a model for in-situ solar wind, 
magnetotail or ground-based magnetometer time series. One reason 
why sLm has not seen wider use here is because it cannot reproduce 
the correlated increments seen for both these types of data and also 
because it models super diffusive (H > 0.5) rather than the observed 
subdiffusive (H < 0.5) behaviour. The term "truncated Levy flight" 
usually indicates standard Levy motion with a finite variance intro- 
duced deliberately by means of a finite range cutoff (c.f. the discussion 
in section 8.4 of Mantegna and Stanley(2000)); however any finite (and 
thus finite-variance) series of sLm must also be naturally truncated, 
but in an uncontrolled fashion (Nakao, 2000). 
Mandelbrot(1995) went on to note that: 

sLm and fBm, however, are far from exhausting the anomalies found 
in nature ... many phenomena exhibit both the Noah and Joseph ef- 
fects and fail to be represented by either sLm or fBm ... One obvious 
bridge, fractional Levy motion, is interesting mathematically, but 
has found no concrete use" . 
Since those words were written, fLm has found applications, notably 
in geophysics (Painter and Patterson, 1994) and telecommunications 
network modelling (Laskin et al., 2002). We here apply it to essentially 
the same need; to compactly describe and unify the "stylized facts" of 
the well-demonstrated Noah and Joseph effects in space plasma physics 
time series (Watkins, 2002). 

3.2. Mathematical definition of fractional Levy motion 

Fractional Levy motion can be defined using a Riemann-Liouville frac- 
tional integral generalising the better-known expression for fractional 
Brownian motion (Voss, 1985). We here adapt the notation of equation 
5 of Laskin et al.(2002), which defines a process W^p: 



draftl3.tex; 2/02/2008; 6:10; p. 4 



Synthesis of scaling exponents: a fractional Levy model 5 



W,At) = ^J2) jf (* - r) {m - l) dW,(r) (1) 

Equation (1) can be unpacked as a summation of Levy stable incre- 
ments dW^(r) each weighted by a response function (t - r)^/ 2 " 1 ). The 
\i parameter describes the power law tail of the pdf of dW which falls off 
as P{x) ~ x~( 1+M "\ 11 = 2 is the special, Gaussian, case corresponding 
to fBm. (3 is the parameter which controls long-range dependence. It 
is well known to be related to the power spectral density S(f) ~ f^ 13 
for fractal processes with finite variance (Voss, 1985), but can also 
be rigorously defined through fractional differentiation in other cases 
(Chechkin and Gonchar, 2000b). 

With ji = 2 and taking in addition (3 = 2 the response function 
becomes unity, giving an uncorrelated random Gaussian walk (WBm). 
Keeping (3 = 2 but allowing /j, to vary in the range to 2 describes 
sLm. fLm is thus in general a process with /x, (3 allowed to vary in the 
range [0 < fi < 2, 1 < (3 < 3] and so forms a bridge between the (3 = 2 
sLm and fi = 2 fBm "axes". fLm thus by construction exhibits both 
the sources of anomalous diffusion identified by Mandelbrot above. 

These limits have corresponding simplified Fractional Kinetic Equa- 
tions (FKE) for the pdf P(W), see section 5.2 of (Zaslavsky, 2002). 
Putting W = W^>(x,t) with (3' = (3/2, WBm is given by the dif- 
fusion equation d\P{W 2 ,i) = 8%(AP(W 2 ,i))i fBm by df P(W 2 ,/?') = 
dl(AP(W 2 ,p>)); and sLm by d\P{W^ 2 ) = d%(AP(W^ 2 )). fLm should 
thus correspond to equation (132) of (Zaslavsky, 2002): 

w p{w ^' ] = ^xY {AP{w ^ ,)] (2) 

All cases have a fixed diffusion constant A. Future work is required to 
establish if this simplified form of equation (127) of (Zaslavsky, 2002), 
the full FKE, can map on to the Fokker-Planck equation of (Hnat 
et al., 2005) or whether the full equation, including fractional drift 
and diffusion terms, is needed. After initial submission of this paper 
we also became aware of the relevance of the work of Milovanov and 
Zelenyi(2001) to the interpretation of fLm as an FKE; see in particular 
their equation (3). 

3.3. Self-similarity, the Hurst exponent and peak scaling 

We now follow Laskin et al.(2002) to show that W^p is indeed an H- 
selfsimilar process. To see this we first put r = cs in (1). We then use 
the fact that the increments dW^{cs) are defined to be l//x self-similar 
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i.e. are equal in distribution ( = ) to c 1 ^dW fl (s). Then 



W^(ct) = c H W^(t) 



(3) 



with a self-similarity parameter H given by 



H = P/2 + l/»-l= [(3/2 - 1/2] + [1/fi] - 1/2 



(4) 



more usually known as the Hurst exponent. Note that we would not 
necessarily expect this equation to hold for more general fractal pro- 
cesses. In the fBm case \x = 2 and for that case only we recover the 
well known expression that (3 = 2H + 1. In the sLm case (3 = 2 and we 
find H = Recently Mandelbrot(2002) has proposed writing 



where he defines a Joseph (long range dependence) exponent J (= 
(3/2 — 1/2) and a Noah (heavy tail) exponent L (=l//x). 

The first property that needs to be shown in a time series for fLm 
to be a candidate model is thus i^-selfsimilarity. This can be tested by 
a number of methods. The first is scaling collapse, which was shown 
for the datasets in our paper by Hnat et al. (2002b; 2003a). 

An fLm model also implies that the pdf of returns i.e P(5X = 0, r) 
will scale with r with exponent also equal to H. This was shown in 
Fig. 2 of (Hnat et al., 2002b). For convenience in figure 1 we show 
a comparison of the scaling regions of the 1 year signals taken from 
the natural time series AE, AU, AL and e. All are seen to scale up to 
approximately 2 6 minutes (~ 1 hour). Caution is however necessary 
because in a natural dataset the moments \5X\ g = \X(t + r) — X(t)\ q 
would be expected to be dominated in the small r limit by the scaling 
of the measurement noise on the differences 5X rather than that of the 
physical variables themselves (Hnat et al., 2005). 

Interestingly, although the exponent needed to rescale the pdfs P(SX, r) 
of differences SX taken from fLm is the "full" extended H = H(p, (3) 
defined in equation (4), the difference pdfs have the same shape they 
would have for an sLm with the same fj, value. This is analogous to the 
way in which fBm retains the same Gaussian distribution as the steps 
from which it is composed, despite their statistical dependence, and is 
why fLm is also known as "linear fractional stable motion" . 



H = J + L-l/2 



(5) 
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Figure 1. Estimation of Hurst exponent H via scaling of peaks P(0) of pdfs of 
differenced time series X(t + r) — X(t) as a function of differencing interval r. Plots 
are for i auroral indices 1978: X = AU(o), X = AL (A) and AE (box) and ii solar 
wind e ({>) for 1995. Plots have been offset vertically for clarity. 



3.4. Structure functions S q and their scaling exponents 
C(q): H AS C(l), WHILE THE PDF of returns gives 
C(-l) = -H 



One may extend the idea of self-similarity expressed by H to the 
generalised q-th order structure functions (Frisch, 1995): 

S q =< \x(t + t) - x(t)\ q > (6) 

where q need not be integer. If a given S q is empirically found to be a 
power law we can then define an exponent from S q ~ T^ q \ 

For a stable self-similar process where all moments are finite = 2, 
i.e. WBm (H = 0.5) or fBm (0 < H < 1), the exponents of the structure 
functions C(q) follow £(q) = qH, as we have checked by simulating 
an fBm using the same fLm algorithm as used for the figures, in the 
fj, = 2 limit. By definition we then have £(1) = H. Additionally in 
these Gaussian (/i = 2) cases Q{2) = 2H, which from Equation (4) then 
implies 2H = (3 — 1. 

The exponent derived from the pdf of returns can be shown to be 
equivalent to £(— 1) [Miriam Forman, private communication, 2002] so 
for self-similar processes (see also our figure 5) the plot of ((q) versus q 
is antisymmetric about q = at least insofar as C( — 1) = — H = — £(1). 
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3.5. Second order moment and J: Pseudo-Gaussian 

BEHAVIOUR OF TRUNCATED LEVY TIME SERIES 

Because of the relation ((q) = qH for WBm and fBm, a complementary 
estimate of the self-similarity parameter H can, for these cases, be 
obtained from from the well-known growth of the standard deviation 
<t(t) of the difference time series 5X(t) with differencing interval r. 
Indeed the growth of a measured a as r 1 / 2 in the case of WBm defines 
diffusive behaviour, a is the square root of variance and thus scales like 
S2, i.e. as t^ -1 )/ 2 , i.e. it follows Mandelbrot's (2002) Joseph exponent 
J (which from (4) will be identical to H in the Gaussian WBm or fBm 
cases) . 

In the case of Levy motion, however, whether ordinary or fractional, 
the qth order moments S q (where q > fi) taken from a set of N 
data points are theoretically infinite as N — ► 00 in contrast to the 
convergence seen for Gaussians. It is thus not a priori obvious how the 
variance of a truncated, finite- N, time series would be expected to scale. 
This is significant because any simulation that we perform of fractional 
Levy motion is effectively one of truncated Levy motion; while a natural 
time series will also have a finite variance in practice. The possible 
relevance of this question to data is clearly illustrated by our Figure (2), 
(see also table 1 of (Hnat et al., 2002b)) in which a for the solar wind 
variable e is seen to scale with an exponent of 0.29 as opposed to the 
values around 0.43 — 0.45 seen for the 3 geomagnetic index quantities. 
Rather than scaling with H, a still appears to be showing pseudo- 
Gaussian behaviour i.e. following J, in that (3 = 1.56 for this time 
series (estimated by wavelet methods) giving J = (1.56 — l)/2 = 0.28. 

The apparent disadvantage of the loss of a second, independent, 
estimate of H seems to be compensated for by the possibility that we 
can use the growth of a to measure (3 i.e. we can effectively use it as a 
measurement of J. On the assumption that a naturally truncated fLm 
describes our data we can build a table (Table I) of the measured (5 
and H values and then predict /x using equation (4). 

On inspecting Table I the first point is that the values of H and 
J are so close in the case of AU that if we assume they are exact 
the predicted [i becomes 2, eliminating fLm as a model for AU. The 
H is sub-diffusive, so fBm would remain a possible candidate model; 
however the observed (Hnat et al., 2003a) difference pdfs P(SX) for AU 
are non-Gaussian, eliminating fBm. The error bars quoted in Table I 
suggest these conclusions may be too harsh. fLm would, however, seem 
more suitable as a model for AE, AL and s. 

As a test we may also consider the values of H and J for solar wind 
B 2 obtained by Hnat et al. (2002a). Their figure 3 gives H = 0.42 in 
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Figure 2. Estimation of exponent J for scaling of the standard deviation a of the 
differenced series versus r for the same quantities as figure 1. Notation as figure 1. 



Table I. Measured values of H (from Fig. 1) and J (from Fig. 2) for natural time 
series, and /j, value predicted from equation (4) on the assumption of naturally 
truncated fLm. All measured values are ±0.02 except J for e which is ±0.03. 



Variable 


Measured H 


Measured J 


Inferred /3 


Inferred L 


Predicted fj, 


AE 


0.47 


0.45 


1.90 


0.52 


1.92 


AU 


0.43 


0.43 


1.86 


0.5 


2 


AL 


0.46 


0.43 


1.86 


0.53 


1.88 


£ 


0.45 


0.29 


1.58 


0.66 


1.51 



our parlance, while they report a scaling exponent for a of 0.28 (i.e. J). 
Inserting this into equation (4) predicts \x = 1.56, which is equivalent 
to the 1/q of their equation (3) (see also their Figure 4) which they 
find to be 1/0.66 = 1.5, encouragingly good agreement. 

3.6. Fractional Levy simulation: Comparison with first 
and second order measures 

We can then now simulate fLm using parameters drawn from natural 
data to see if the inferences we have drawn above are indeed consistent, 
and to qualify fLm as at least a possible proxy for these time series. We 
use the published algorithm of (Wu et al., 2004). This has the advantage 
of being linked more closely to the definition of fLm from equation 
(1) than the (faster) approach of replacing (Chechkin and Gonchar, 
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2000b) a Gaussian random number generator by a Levy generator in 
otherwise standard Fourier filter methods (Voss, 1985). A comparison 
of these two approaches will be reported in a future paper. 

We show simulation results for synthetic AL and e time series. These 
were specified by the ordered pairs ((3,fi) of (1.86,1.88) and (1.58,1.51) 
respectively. The P{5X = 0, r) scaling for both series (Figure 3) is 
seen to follow H as we expect, so both model series have very similar 
measured H values, as we also saw in their natural counterparts (Figure 
(1)). Conversely, for finite samples of fLm, however, modelling AL and 
e we see from Figure (4) that rather than following r 1///1 (= t l ), the a 
measured on the difference time series 5X still grows as r^ _1 )/ 2 (= t j ) 
i.e. it does, as postulated in subsection 3.5, measure J rather than L. 

This effect requires some discussion. It seems to be a further man- 
ifestation of the "pseudo-Gaussian" behaviour of truncated standard 
Levy motion (Chechkin and Gonchar, 2000a), and known (Nakao, 2000) 
to be responsible for the result ((2) = 1 in that case (see also Fig- 
ure 5). Our simulations have clearly demonstrated that it generalises 
to the long-range dependent fLm case i.e. that in general for fLm 
C(2)/2 = J = ((3 — l)/2. This conclusion is most clearly supported 
by Figure 5 where the ((2) value can be read off as following this 
relation over the range [3 = 1.5 to 2.5. The agreement is poorer at 
smaller (3 values tested. We currently think this reflects known difficul- 
ties with accurately simulating strongly anticorrelated fLm (Chechkin 
and Gonchar, 2000b). The effect has previously been remarked on in the 
truncated standard Levy paradigm; for example the S&P 500 financial 
time series, depicted by (Mantegna and Stanley, 2000) where [3 = 2 
(their Fig. 11. 4. a) so a grows as t 1 / 2 (their Fig. 11.3a), in contrast to 
an H value from peak scaling of 0.71 (their Fig. 9.3). 

In the multifractal modelling community the power spectrum has 
long been seen as only just one of several ways of measuring £(2). For 
this reason a difference in the value of £(1) ^ ((2)/2 has sometimes been 
claimed as direct evidence of the inapplicability of any additive model 
and thus the immediate need for a multiplicative model (Schertzer and 
Lovejoy, 1987). 

Conversely our result would seem to suggest that any truncated 
stable additive model other than the fBm/WBm limiting cases is likely 
to show C(l) + C(2)/2, and C(2)/2 = J = ((3 — l)/2 without the need 
for a multiplicative model. This may be understood as being because 
truncated Levy motion, whether standard or fractional, behaves as a 
bifractal (Nakao, 2000). There may be natural time series where addi- 
tive fLm is actually the most natural model, or at least an economical 
and easily specified one. 
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Figure 3. Estimation of H via scaling of peaks P(0) of pdfs of differenced model 
fractional Levy motion time series X(t + t) — X(t) as a function of differencing 
interval r. Plots are for i) a synthetic AL(o) time series and ii) a series of synthetic 
solar wind X — e (<>)• Plots have been offset vertically for clarity. 



3.7. C PLOTS AND THE MULTIFRACTALITY OF TRUNCATED LEVY 
MOTIONS 

At this point it may be objected that we have not tested any predictions 
of the fLm model against the behaviour of natural time series other 
than those properties used to specify it. Our first additional check is 
thus to examine the multi-affine behaviour seen in the data and the 
model using the % plots" defined in section 3.4. Such a plot, showing 
scaling exponent Q{q) versus moment q is shown for the data in figure 
(6). Interestingly AU most resembles a "classic" multifractal, in that 
the points £(17) lie on a curve rather than a straight or broken line 
(Frisch, 1995). However AE, or at least AL, have ( which arguably 
flattens out near 1 for higher moments, e intriguingly even seems to 
fall as m increases. This behavior is qualitatively similar to that seen 
for our simulated AL and e time series, whose Q{q) plots are superposed 
on the figure. In particular a change in the range of r over which the 
simulated AL structure functions are taken to be power laws is enough 
to encompass the observed range of £ plots for e. More detailed com- 
parison is at present prevented by the difficulty of obtaining accurate 
values of S q for high moments-an issue also afflicting analysis of real 
data. 
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Figure 4- Estimation of J by scaling of the standard deviation a of the differenced 
simulated (and thus truncated) series versus r for the same quantities as Figure (2). 
Notation as in (2). 



A significant body of data and models now exists for the problem 
of solar wind and magnetic index scaling. We have here suggested 
a complementary approach, motivated in particular by the need to 
i) reconcile differing estimates of scaling exponents (in hindsight the 
Joseph and Hurst exponents J and H); ii) model subdiffusive behaviour 
(H < 0.5); and iii) model long-ranged correlation ([3 / 2). We proposed 
the use of a simple and economical model: fractional Levy motion, to 
describe the scaling of the above quantities. Initial consistency checks 
with respect to the distribution of returns and the scaling of standard 
deviation support the use of fLm, and the multi-affine "zeta plots" 
are more qualitatively similar. Importantly we find that the degree of 
similarity between model solar wind e and the model AL index does 
indeed depend on the moment order at which comparison is made, 
but that this does not, however, require a multiplicative process to 
explain it. The difference can, rather, be understood as coming from 
the bifractality of a truncated fractional Levy motion. This explains 
why some measures such as H from the distribution of returns or pdf 
rescaling are much closer to each other than, for example, the <r-based 
exponent (which we found to measure J, not H). 



4. Conclusions 



draftl3.tex; 2/02/2008; 6:10; p. 12 



Synthesis of scaling exponents: a fractional Levy model 



13 




_1 1 ■ ' • • 

—10 12 3 4 

q 



Figure 5. Zeta plots for simulated (and thus truncated) fLm with /i fixed at 1.5 and 
(3 ranging from 1.1 to 2.5. The relation £(2) = (3 — 1 is seen to be well satisfied for 
/3 > 1.5. 

The present paper has been mainly concerned with the modelling 
of measured quantities rather than the extent to which they are artifi- 
cial. For geomagnetic indices and other constructed quantitities like e, 
however, the extent to which scaling behaviour could be an artefact of 
the construction method is an important issue. We are aware of some 
progress in studying this problem (e.g. Edwards (2001), Weigel and 
Baker (2003)), more will be needed. Further work is also underway to 
test the predictions of the fLm model against other scaling studies such 
as the cited burst lifetime and spreading exponent investigations. 
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